Quantum fluctuations in Larkin-Ovchinnikov-Fulde-Ferrell superconductors 



K. V. Samokhin 1 and M. S. Mar'enko 1 ' 2 
1 Department of Physics, Brock University, St. Catharines, Ontario, Canada L2S 3A1 
2 Department of Physics and Astronomy, Hofstra University, Hempstead, New York 11,54-9, USA 

(Dated: February 2, 2008) 

We study the superconducting order parameter fluctuations near the phase transition into the 
Larkin-Ovchinnikov-Fulde-Ferrell state in the clean limit at zero temperature. In contrast to the 
usual normal metal-to-uniform superconductor phase transition, the fluctuation corrections are dom- 
inated by the modes with the wave vectors away from the origin. We find that the superconducting 
fluctuations lead to a divergent spin susceptibility and a breakdown of the Fermi-liquid behavior at 
the quantum critical point. 
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I. INTRODUCTION 



Magnetic field suppresses superconductivity, regardless 
of the pairing symmetry, via the coupling to the orbital 
motion of electrons in the Cooper pairsi In spin-singlet 
superconductors, the pairs are also broken by the Zee- 
man interaction of electron spins with an applied field 
H, or by the exchange interaction with localized spins 
in a magnetic crystal, which is known as the paramag- 
netic, or Pauli, mechanism. If the orbital effects are ne- 
glected, then as shown by Larkin and Ovchinnikov^ and 
Fulde and Ferrcll 3 (LOFF), the competition between the 
paramagnetic pair breaking and the condensation energy 
results in the formation at low temperatures of a peculiar 
non-uniform superconducting state with a periodic mod- 
ulation of the order parameter, whose critical field H c 
exceeds the Clogston-Chandrasekhar paramagnetic limit 
for a uniform state4*£> The superconducting order pa- 
rameter in the simplest realizations of the LOFF state is 
either a single plane wave: rj(r) = ri e lqcr , or a super- 
position of two plane waves: rj(r) = 770 cosq c r, where q c 
is the wave vector of the LOFF instability. In general, 
the order parameter structure can be more complicated 
and is determined by minimizing the nonlinear Ginzburg- 
Landau free energy. 

For a long time the LOFF state had been considered 
a theoretical curiosity, because its experimental detec- 
tion required the fulfillment of some rather stringent 
conditions. First, the orbital effects are detrimental to 
the LOFF state and therefore should be weak enough. 
The relative importance of the orbital and spin pair- 
breaking mechanisms is measured by the Maki param- 
eter olm — \Z2Hjg /Hccfi where i2^2 i s *he upper crit- 
ical field in the absence of spin interactions, and Hqc 
is the Clogston-Chandrasekhar critical field. In the pure 
paramagnetic limit, the orbital pair-breaking is absent al- 
together and a a/ = 00. The orbital effects were included 
in the LOFF model with a spherical Fermi surface by 
Gruenberg and Guntherji who found that at T = and 
a M ^ 1-8 the order parameter below the upper critical 
field is modulated along the direction of the applied field. 
The coordinate dependence of the pair wave function in 



the transverse directions is described by the lowest Lan- 
dau level, i.e. is the same as in the pure orbital casei^ 
There is another phase transition at a lower field, either 
into the usual mixed state or into the uniform supercon- 
ducting state at gym = 00, resulting in the appearance 
of a characteristic wedge-like region in the H — T phase 
diagram at low temperatures and high fields. In most 
"classical" superconductors, however, aM < 1, so the 
orbital pair-breaking dominates and the LOFF state is 
never realized. 

One possible way to reduce the orbital effects was pro- 
posed by Bulaevskiijifl who pointed out that in a layered 
superconductor with the electron orbital motion confined 
to the layers, the Maki parameter depends on the angle 
9 between the direction of H and the layers, making the 
paramagnetic effects dominant in a narrow angle interval 
near 9 = 0. As 9 approaches zero, the system undergoes 
a series of phase transitions between the LOFF states 
corresponding to successive higher Landau levels. The 
full H — T phase diagram was worked out in Ref . [H At 
the parallel field orientation there is no orbital effects, all 
pair breaking is entirely paramagnetic, and the region of 
existence of the non-uniform state turns out to be larger 
than in the isotropic 3D case, see also Ref. [l2|- The same 
ideas also apply to thin superconducting films, 13 or to 
surface superconductors^, in parallel fields. 

Another obstacle to the experimental realization of the 
non-uniform state is its sensitivity to the presence of dis- 
order. It was found by Aslamazo^ii^ in the isotropic case 
that the LOFF critical field decreases rapidly with in- 
creasing non-magnetic impurity scattering and eventu- 
ally becomes smaller than Hcc: resulting in the restora- 
tion of a first-order phase transition into the uniform su- 
perconducting state. Later the analysis was extended to 
the layered case in Ref. Il6l with essentially the same 
conclusions. 

Thus the LOFF state can potentially be observed 
only if the superconductor is both paramagnetically lim- 
ited and sufficiently clean. These requirements can 
be simultaneously met in heavy-fermion compounds, in 

which H^) is inherently high due to a short coherence 
length. Earlier candidates for hosting the LOFF state in- 
cluded UPd 2 Al 3f iLi£ UBei 3r i2 and CeRu 2 ^ The odds 
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of finding the LOFF state are even greater in quasi- 
low-dimensional superconductors, in which the orbital 
pair breaking is reducediSi'SI* 2 ^ 2 ^ Several experiments 
on organioSiS 5 ^ and cuprate^ superconductors have re- 
vealed the features in the the H — T phase diagram, such 
as an upturn of the upper critical field and the presence 
of an additional phase transition below H C 2 at low tem- 
peratures, that could be interpreted as signatures of the 
LOFF state. Similar features have been recently reported 
in another heavy- fermion compound, CeCoIns ^2*22 In all 
the cases mentioned above, the spin splitting of the elec- 
tron energies in the Cooper pairs was due to the Zccman 
interaction with an applied magnetic field. It was ar- 
gued in Rcf. that the LOFF state can be created by 
the intrinsic exchange band splitting in the ferromagnetic 
superconductor RuS^GdC^Os- 

Another intriguing possibility of the experimental re- 
alization of the LOFF state has been discussed very re- 
cently in the context of ultracold atomic Fermi gases, 
such as 40 K and 6 Li. By making the populations of atoms 
in two different hyperfine states unequal, one controls 
the mismatch between their Fermi surfaces When 
the pairing interaction between the two fermion species 
is turned on, the system becomes formally equivalent 
to a neutral superconductor in a Zeeman field. Due to 
the absence of both the orbital effects and impurities, 
this seems to be the most promising setup to study the 
paramagnetic pair breaking, including the non-uniform 
states!??'? 4 !?^?*?!?! 

Despite the lack of unambiguous experimental evi- 
dence, the LOFF state has remained a subject of inten- 
sive theoretical investigations in the past decades. In 
addition to the studie s cited above, we would like to 
mention Refs. in which the Ginzburg-Landau 

theory was developed in the pure paramagnetic case in 
the vicinity of the tricritical point in the H — T phase 
diagram, see Fig. 2 below, where the sign change of 
the second-order gradient term in the free energy sig- 
nals the onset of the non-uniform instability. The vor- 
tex st ructure in t he mixed LOFF state was studied in 
Refs. l4(l4ll42l48L The LOFF model has also been ex- 
tended to unconventional, in particular c£-wave, pairing 
symmetries ^^MJ^ An analysis of the spatial structure 
of the LOFF state immediately below the upper critical 
field in the isotropic 3D case was done by Larkin and 
Ovchinnikovj 2 . who showed that it is the "striped" phase 
with r](r) = rjo cos q c r that is energetically favored at 
T = 0. The zero-temperature phase transition from the 
normal state was found to be second order, but becomes 
first order as temperature increases l iL224 9 . On the other 
hand, it was argued in Refs. I50ll5ll that the phase transi- 
tion is always first order below the tricritical point, and 
that the order parameter at T = is represented by a 
sum of three cosines. 

Another open question concerns the nature of the lower 
phase transition separating the LOFF state from the 
conventional uniform superconducting state. The only 
cases studied so far assumed a one-dimensional periodic- 



ity of the order parameter in a purely paramagnetic and 
isotropic 2D2fl* 4 ii or 3D3& system. As the field decreases, 
the non-linear effects add higher harmonics to the LOFF 
state, which starts to resemble a periodic array of Bloch 
domain walls separating the regions where the order pa- 
rameter is almost uniform. When one approaches the 
lower critical field, the period of the domain structure 
diverges, indicating a second order phase transition. 

While an extensive literature exists about the mean- 
field properties of the LOFF state, the superconducting 
fluctuation effects have received comparatively little at- 
tention. The phase fluctuations of the non-uniform order 
parameter were considered by Shimahara, 52 who found 
that they are able to destroy even the quasi-long-range 
order for the striped LOFF states in the isotropic 2D case 
at T > 0. He also conjectured that in the isotropic 3D 
case the long-range order is replaced at finite tempera- 
tures by a quasi-long-range order. This is consistent with 
the findings of Ref. 0^ where it was shown that the ther- 
mal fluctuations suppress the second order phase transi- 
tion into the LOFF state in spatially isotropic systems. 
These effects are analogous to the fluctuation-driven de- 
struction of crystalline order with one-dimensional den- 
sity modulation^ In general, one can expect that, since 
the wave vectors of important fluctuating modes in the 
isotropic LOFF state are close to a sphere (or a circle in 
2D) of radius \q c \ ^ 0, the fluctuation effects on observ- 
able quantities will be considerably enhanced compared 
to the uniform case due to the increased phase volume 
of the fluctuations^ The fluctuation effects might still 
be significant even when the degeneracy manifold of the 
LOFF states is reduced to a set of isolated points in the 
momentum space: It was recently argued in Ref. l5fil that 
the thermal fluctuations in quasi-2D d-wave supercon- 
ductors are strong enough for the LOFF phase transition 
to become of first order. 

In all the works mentioned above only finite tempera- 
tures were considered, in which case the order parameter 
fluctuations are predominantly classical. Formally, the 
classical limit corresponds to setting the frequency in the 
fluctuation propagator to zero, see Sec. HT1 below. The 
focus of the present work is on the fluctuation effects at 
T = above the quantum phase transition from the nor- 
mal state to the LOFF state driven by an external mag- 
netic field. In this case, the dynamic nature of the super- 
conducting fluctuations cannot be neglected. We assume 
that the system can be described by the Bardeen-Cooper- 
Schrieffcr (BCS) model and also that the quantum LOFF 
transition is of second order. We do not include impu- 
rities and the orbital effects, expecting our results to be 
applicable either to paramagnetically limited supercon- 
ductors, or to the Fermi gases of ultracold atoms. We 
would like to note that while thermal fluctuations in su- 
perconductors have been actively studied for a long time, 
see Ref. 15^. the quantum fluctuations at low tempera- 
tures have only recently become a subject of theoretical 
investigation. 58 ! 59 ! 60 ! 61 ! 62 

The paper is organized as follows: In Sec. [B] we de- 



rive the general expression for the fluctuation propagator 
in the normal state above the LOFF phase transition. 
We consider both the isotropic case, in which the LOFF 
states are infinitely degenerate in the momentum space, 
and the generic case, in which the degeneracy is lifted due 
to the band structure and/or the gap anisotropy. In Sec. 
IIIII we calculate the quantum fluctuation corrections to 
the spin susceptibility and to the decay rate of fermionic 
quasiparticles. 



II. LOFF FLUCTUATION PROPAGATOR 



3 




FIG. 1: The diagrammatic representation of the fluctuation 
propagator. The solid lines are the Green's functions of 
fermions in the normal state, the bold dots denote the BCS 
pairing potential. 



We consider a clean spin-singlet BCS superconductor 
in an external magnetic field H . The coupling of the 
electron charges to the vector potential is neglected, so 
that the superconductivity is affected by the field only 
through the Zceman splitting of the single-particle bands. 
The Hamiltonian is given by 

H = ^(6=<W ~ h(J 3ya f3)cl a C k p 
k 

+ F W9) c fe 1+ f,T c lfe 1+ i4 c -fe 2 +fd c fe2+§,t' (!) 
q,fci,2 

The first term here is the free-fermion part, where = 
ffc — fx, £fc is the band dispersion, fi is the chemical poten- 
tial, a, (3 =1, 1 is the spin projection on the quantization 
axis along if, h = hbH is the Zeeman field, fiB is the 
Bohr magneton, and 03 is the Pauli matrix (we use the 
units in which K = ks = L and assume that the Landc 
factor g = 2). The Hamiltonian (JIJ can also be applied 
to a ferromagnetic superconductor in zero applied field, 
in which case the electron bands are split due to the ex- 
change interaction with the magnetically ordered local- 
ized spins. 

The second term in Eq. is the pairing interaction, 
which is effective only in the vicinity of the Fermi surface 
defined by the equation = 0, i.e. at |£±fc 12 +q/2| < 
u max , where w max is the BCS energy cutoff. We make a 
simplifying assumption that the interaction matrix can 
be factorized: 

Vfcifc2 (<?) = -Hq)^ 0k 2 , (2) 

where cf>k = 4>~k is the symmetry factor, which is nonzero 
only inside the BCS shell, i.e. at |£fc| < w max . The 
symmetry factor is assumed to be real and normalized: 
((jyj.) = 1, where the angular brackets stand for the Fermi- 
surface average. In the group-theoretical language, 4>k is 
the basis function of an even one-dimensional irreducible 
representation T of the point group Q of the crystal, 
which can have zeros, symmetry-imposed or accidental, 
somewhere on the Fermi surface. The pairing is said to 
be conventional if T is the unity representation, and un- 
conventional otherwise & To make sure that the energies 
of all four fermions participating in the BCS interaction 
are less than <x> max , one has to further assume that the 



function A(q) is nonzero only if \q\ < g max ~ tL> max /vp, 
where vp is the Fermi velocity. In the calculations below 
we replace X(q) by a coupling constant A > 0, introduc- 
ing an explicit momentum cutoff in q-integrals if needed. 

The order parameter dynamics in the normal state is 
described by the fluctuation propagator 57 



where v m — 2-kttiT is the bosonic Matsubara frequency, 
and C{q,v m ) is the particle-particle propagator (the 
Cooperon), see Fig. 1. Calculating the diagrams, we 
obtain: 

^£-W> = (2) 

+ (*H..(i + ^)), (4) 

where Np is the density of states per one spin projec- 
tion at the Fermi level, ^(x) is the digamma function, 
*(l/2) = -ln(4e c ), C ~ 0.577 is Euler's constant, 
T c o = (2e c /7r)(jj max e -1 / ArFA is the zero-field critical tem- 
perature of the uniform superconducting state, 

W k = £ fe+§ - £ fe -s - 2h = v k q -2h + 0(q 3 ), (5) 

and Vf, = Vfe^fc is the quasiparticle velocity at the Fermi 
surface. The cutoff w max has been eliminated by adding 
and subtracting the Cooperon at h = v m = q = 0. 

The expression (0} is valid at all temperatures, for 
any pairing symmetry and band structure [we will keep 
only the linear in q term in the expansion JSJ, assum- 
ing that the band structure is such that the higher-order 
terms are negligible]. The solution T(q,h) of the equa- 
tion £ _1 (q, 0) = determines the temperature at which 
the superconducting instability with the wave vector q 
develops in a given field h. Setting v m = q = T = 0, 
one finds that the second-order quantum phase tran- 
sition into a uniform superconducting state occurs at 
h = (n/2e c )T c0 = A /2 ~ 0.88T c0 , where A is the 
BCS gap at T = 0. In general, the critical temperature 
vs field T c (/i), or inversely the critical field vs temper- 
ature h c (T), can be found by maximizing T(q,h) with 
respect to q. According to Refs. 00, in a clean isotropic 
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FIG. 2: The generic phase diagram of the LOFF supercon- 
ductor. The LOFF state appears at H = H c (solid line) at 
temperatures below the tricritical point (T* , H*). The dashed 
lines correspond to the normal metal-to-uniform superconduc- 
tor phase transitions: the first-order Clogston-Chandrasekhar 
transition at Hcc (dashed line), and the second-order transi- 
tion at Ho (long-dashed line). The critical field of the LOFF 
state-to-uniform superconductor transition is not shown. 



superconductor at T < T* ~ 0.56T c o the maximum of 
the critical field is achieved at q c ^ 0. The generic phase 
diagram of a LOFF superconductor is sketched in Fig. 2. 

In the vicinity of the critical field h c (T), the most di- 
vergent contributions to physical quantities come from 
the low-frequency fluctuations, so the inverse fluctuation 
propagator can be expanded in powers of v m : 

-^C-\q, v m ) = A(q, h) + A(q, h)^ + 0(v 2 J, (6) 



where 



A(q,h) 



A(q,h) 



4 Re* 



2ttT 



1 iW k \ 

2 4ttTJ 
T 



In- 



T, 



c0 



iWk \ 
AttTJ 



(7) 

(8) 



We focus on the fluctuation effects at T = 0, when 
the explicit temperature dependence of the expressions 
O and (JSJ) can be eliminated by using the asymptotic 
form of the digamma function: ~^>(x) = In a; + 0(x^ 1 ) at 
x — ► oo (Ref. l64j) : 



A(q,h) = In ^- +F(Q), 
A(q,h) = F(Q), 
where Q = q/2h, and 

F(Q) = (<fthi\v H Q-l\), 

HQ) = «(tf.S(w.Q-i)) 



(9) 
(10) 

(11) 

(12) 



If the function F(Q) has a minimum at Q — Q c , then 
the upper critical field is given by h c — hoe~ p ^ c ' , and 
the equilibrium wave vector of the LOFF structure is 
q c = 2h c Q c . 

In the standard theory of superconducting fluctua- 
tions, see Ref. l57l it is assumed that the maximum of the 
critical temperature, or of the critical field, is achieved for 
the uniform superconducting state, which makes it pos- 
sible to expand £ _1 (q, v rn ) in the vicinity of the origin in 
the momentum space. In contrast, the most important 
critical fluctuations in the LOFF state have the wave vec- 
tors near q c ^ 0. Assuming that F(Q) can be expanded 
in the Taylor series near the minimum, we obtain the 
following expression for the fluctuation propagator: 



£{q,v m ) 
1 



where 



N F e + j\v m \ + K lj (q l - q c ,i)(qj - q c ,j) ' 



h — h r 



(13) 



(14) 



measures the distance to the quantum critical point, 



HQc) 

2h c 



(15) 



In should 



and Kn = (l/8hl)ViV s F(Q c ), i,j = x,y,x. 
be noted that in some cases the frequency expansion © 
of the inverse fluctuation propagator does not exist, see 
Sec. UTTl below. 



A. Isotropic 3D case 

Explicit expressions for 7 and K%j can only be ob- 
tained in few cases, including a 3D parabolic band, 
£fc = k 2 /2m — /i, with isotropic pairing. In this case 
a straightforward integration in Eqs. (|11I12|I gives 



HQ) 
HQ) 



1 



I11V 



1 

2x 



In 



1 



2.i; 



9(x-l) 



1, (16) 
(17) 



where x = vpQ = Vpq/2h. The function (|16f) has a min- 
imum at x = x c ~ 1.20, with F(Q C ) ~ —0.41, see Fig. 
3. Thus the quantum phase transition occurs at h c ~ 
1.51/io — 0.75Ao, and at h < h c the superconducting 
order parameter is spatially modulated, with the wave 
vector q c — 2x c h c /vF — 0.51^ , where £0 = vf/2ttT c q is 
the BCS coherence length. The LOFF critical field h c ex- 
ceeds not only ho, but also the Clogston-Chandrasekhar 
field hcc = which corresponds to a first-order 

phase transition between the normal and the uniform su- 
perconducting states*^ see Fig. 2. 
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FIG. 3: The momentum dependence of £~ 1 (q,0) in the 
isotropic 3D case, Eq. 1161 (solid line), and the isotropic 
2D case, Eq. IJ220 (dashed line); x = vp\q\/2h. 



In the isotropic case the critical field of the LOFF state 
does not depend on the direction of q, and the quantum- 
critical fluctuation propagator i|13|) takes the following 
form: 



r( s_J 1 

V ~ N F e + j\u m \ + K(\q\ — q c ) 2 ' 



(18) 



where 7 ~ 0.65/h c and K ~ 0.28v 2 F /h 2 c . The static 
limit of Eq. (|18|l has the same form as the propaga- 
tor of classical fluctuations of the order parameter asso- 
ciated with the crystallization transition in an isotropic 
liquid.— Note also that the minimum of C^ 1 (q,0) re- 
mains infinitely degenerate even for an ellipsoidal Fermi 
surface, = J2i kf/2mi — k F /2m, where rrn = m//if are 
the effective masses, i = x,y, z. In this case, a change of 
variables reduces the fluctuation propagator to the form 

l(Tg)l. in which \q\ - 



Therefore, 



the critical field has the same value as in the isotropic 
case, the only difference being that the degeneracy man- 
ifold in the momentum space is now ellipsoidal rather 
than spherical. 



parameter just below the critical field can be represented 
as a linear combination of the plane waves: 



r)(r) 



N q 
a=l 



(19) 



The complex coefficients ry a , which determine the spa- 
tial structure of the LOFF phase, are found by mini- 
mizing the Ginzburg-Landau free energy. If the minima 
of A(q, h) are well-separated then the fluctuation modes 
near different qi°^ can be treated independently. 

Let us consider for concreteness a tetragonal crystal 
with Q = Di/j, in which case there can be as many as 
sixteen degenerate minima of A(q,h). This number is 

severely reduced if one assumes that q^ are along the 
highest symmetry directions. For a 3D band this means 



that N n 



- A, with qc 1 ' 2 ^ — ±q c z. Near the minima, one 
can write A lj2 (q, h) = e + K\\(q z =F q c ) 2 + K ± {ql + q 2 y ). 
On the other hand, if there are 2D bands in the system, 
then the lowest possible number of the minima is four, 



located for instance at q^c^ = ±q c %) Qc'"*' = ±<?ej/- 
Then, A 1 (q,h) = e + K\\ (q x - q c ) 2 + K ± q 2 , etc. Without 
any loss of generality, we can assume that Kn = K± = K . 
so that the fluctuation propagator near the ath minimum 
can be written in the following form: 



(2,4) 



£ a (<7, V m ) 



1 



1 



Nf , I I . 7v- / (a) 

e + 7|z/ m |+A [q-qh 



(20) 



This expression, in which 7 and K should be treated as 
phenomenological constants, is applicable in the generic 
case of a 3D or 2D crystalline superconductor with arbi- 
trary band structure and pairing symmetry. 

We would like to note the formal similarity between our 
superconducting fluctuation propagators (|18fl and l|20() 
and the propagators of magnetic fluctuations in itinerant 
helical ferromagnets— and high-T c superconductors, 66 
respectively. 



B. Generic band structure 

The infinite degeneracy of the LOFF states is an arti- 
fact of the parabolic band approximation. It will be lifted 
in the case of a general band structure or an anisotropic 
gap symmetry. Since A(g~ 1 q, h) — A(q,h), where g is 
an arbitrary element of the point group G, the minima 
of the inverse fluctuation propagator in the momentum 
space form a "star", i.e. a set of N q isolated points, 

{qi^}, which is invariant under all operations from Q. 
Assuming that inversion is present in the point group, 
N q can be as low as two and as high as the total num- 
ber of the group elements in Q. The equilibrium order 



C. Isotropic 2D case 

The fluctuation propagator does not have the simple 
form l|13f) if 7 is either zero or infinity. According to Eq. 
(|12fl . the former possibility occurs if the surface VkQ c = 1 
does not intersect the Fermi surface, or if it does then 
4>k accidentally vanishes on the intersection line. In ei- 
ther case one would have to go to higher orders of the 
frequency expansion. We have not been able to find an 
explicit example of the band structure for which this hap- 
pens. 

If 7 = 00 then the Taylor expansion in powers of | v m \ 
fails, and one to take the low-temperature limit directly 
in Eq. (J3J. Assuming as before that the temperature is 



6 



the smallest energy scale in the system, one obtains: 



III. FLUCTUATION CORRECTIONS 



N F h 



2h 



(21) 



In order to recover the expressions (|9I1U|) from this, one 
has to replace \v m \ — » \v m \ + + . 

One can check that the frequency expansion fails in the 
case of the isotropic 2D band with = (k 2 + k 2 )/2m — /i 
and (f>k = 1. For the order parameter modulated in the 
xy plane we have 



A. Free energy and spin susceptibility 

If the quantum LOFF transition in the isotropic 3D 
case is first order jSLSi then our theory is not applicable. 
However, even in that case it is still instructive to do 
the calculations using the fluctuation propagator (|TB)l in 
order to highlight the differences with the generic case. 
The fluctuation correction to the free energy in the nor- 
mal state at T = is given by 



dv 



SF = 2 Y,L i ln/r1 ^-) 



(25) 



F(Q) = Re In 

HQ) = 1 



1 



VT^x 2 



1 



2 

(x-1) 



(22) 
(23) 



where x — v F Q — v F q/2h. The function F(Q) has a 
non-analytical minimum at x = 1, with F(Q C ) = — In 2, 
see Fig. 3. Therefore the zero-temperature critical field 
is h c = 2ho, and the LOFF wave vector is q c = Aho/v F . 
Since F diverges at the critical point, one has to use Eq. 
B21J1 . with the following result: 



1 



h 



—C (q, v m ) = In — + T 



( v F \q\ 



h 



V 2h ' 2h 



(24) 



where 



T(x, y) = Re In 



1 + iy + y/{l + iy)< 



We see that the retarded fluctuation propagator C R (q, v) 
obtained from Eq. (|24l) has a branch cut instead of a 
simple pole. 

The non-analyticity of the inverse fluctuation propa- 
gator persists even if the Fermi surface is a corrugated 
cylinder. For the quasi-2D band described by = 
(k 2 + k 2 )/2m — t cos k z d— \x (t <S fi), with isotropic pair- 
ing, one can show that the deepest minimum of A(q, h) 



is achieved for q 



12 



Then, F(Q) and F(Q) are given 



by the same expressions (122t and (|23|l as in the 2D case, 
but with x — td\q z \/2h. The divergence of F at x = 1 
again signals the failure of the expansion of C^ 1 {q,v m ) 
in powers of frequency. 

Below we neglect these complications and assume that 
the fluctuation propagator has either the form 118|) or 
l|2*U)l. This does not seem to be very restrictive, espe- 
cially since a well-defined frequency expansion and ana- 
lytical momentum dependence will be restored in realis- 
tic layered superconductors by the Fermi-surface or gap 
anisotropy, or by disorder. 



The momentum integration is restricted to |q| < g ma x, 
see Sec. [n] In addition, the ultraviolet cutoff i/ max ~ 
w maj[ is introduced to guarantee the convergence of the 
frequency integral. This cutoff can be extended to in- 
finity when calculating the correction to the magnetic 
susceptibility: 

d 2 1 d 2 

S ^ = -8IP SF = - N -HlO? SR 

In the isotropic 3D case, using the fluctuation propa- 
gator H18[) . we have at e — > 0: 



N F f d 3 q 1 
KlH 2 J (2ir) s e + K(\q\-q c 
N F ql 1 



2ir 2 ~/H 2 y/K yft 



(26) 



(the main contribution to the integral comes from the 
vicinity of q c <C g max ). To estimate the magnitude of the 
correction, we compare it to the Pauli spin susceptibility 
in the normal state, xp = 2ji 2 B Np: 



^ . 1.21 



Xp 



(-F 



H-H r 



1/2 



(27) 



where e F = kp/2m is the Fermi energy. Although this 
expression is divergent at the quantum critical point, 
the size of the fluctuation correction at any nonzero e 
is small because of the factor (A /e F ) 2 . The width of 
the fluctuation region in this case can be estimated as 
(H - H c )/H c - (A /e F ) 4 . Note also that the field- 
dependent fluctuation contribution to the magnetization, 



SM 



d_ 



SF 



1 d 

-Np — -5F, 
H c de 



is of the type SM oc y/e, and is not singular at e — > 0. 
In the generic 3D case, using Eq. (1201) we obtain: 



Np 

^H 2 ^ 



d 3 q 



(2tt) 2 



K q-q, 



(a) 



(28) 
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Then, for each minimum q c a ' , the integral can be calcu- 
lated by shifting the integration variable, q — q c a ' = p: 



N F \ - 



d 3 p 



1 



N F q n 



2ir 3 H 2 jK 



(2tt) 3 e + Ap 2 
1 - 



2q„ 




A", 



(29) 



which is not singular at e — > 0, so is the correction to the 
magnetization. 

In contrast, in the generic 2D case with isolated min- 
ima the reduced dimensionality of the momentum inte- 
gral leads to a logarithmic singularity in Sx- 



$X 



E 



d 2 P 



i 



(2tt) 2 e + Kp 2 



N F 



Att 2 H 2 ^K 



ln ^gmax 



(30) 



Using 7 ~ l//i c and K ~ Vp/h 2 , we have the following 
estimate for the quantum fluctuation correction to the 
susceptibility: 



fry 



A . H c 
In . 

e F H-H c 



(31) 



Our analysis can be easily extended to the case of the 
normal metal-to-uniform superconductor transition (one 
should remember though that in a paramagnetically- 
limited clean superconductor this transition does not ex- 
ist, since the LOFF instability always preempts the uni- 
form instability at T — 0). Formally setting q c = in the 
fluctuation propagator, we see that the infinitely degen- 
erate case (|26|) is never realized and the correction to \ 
in three dimensions is non-singular, see Eq. (|29() . In the 
2D case, one would have the logarithmically divergent 
correction (|30|l . 



B. Quasiparticle decay rate 

The lowest-order contribution to the self-energy of 
spin-up fermions due to the superconducting fluctuations 
in the normal state is given by 

E T (fc,w„) = -T^^£(<7,0 

rn q 

xG x (-k + q,-Lj n + v m ), (32) 

see Fig. 4, where Gi(k,u> n ) — (iu> n — — /i) 1 1S the 
Green's function of spin-down fermions (the pairing is 
assumed to be isotropic). At T — we obtain for the 
quasiparticle decay rate at the spin-up Fermi surface, i.e. 
for k satisfying = h: 




*,co n ,t 



FIG. 4: The contribution of superconducting fluctuations to 
the fermion self-energy, Eq. The solid lines are the 

single-particle Green's functions. The wavy line is the fluctu- 
ation propagator C(q,v m )- 



where k is the direction of the Fermi momentum, 
£ R (q, v) is the retarded fluctuation propagator, and the 
integration is restricted to the region in the q-space de- 
fined by the condition < Wk < u>. Below we calculate 
the decay rate in the limit uj — > at the quantum critical 
point, i.e. at e = 0. 

In the isotropic 3D case the fluctuation propagator is 
given by Eq. I|18|) . It is convenient to choose the polar 
axis in the q-space along = v F k (we neglect the dif- 
ference between the quasiparticle velocity at the spin-up 
Fermi surface and the Fermi velocity v F , which is of the 
order of h c /e F ). It is easy to see that the decay rate in 
this case does not depend on k: 



7 



4ir 2 N F J_! 
q 2 dq 



K 2 (q-q c r + 1 2 (LU-W k ) 2 



(s = cos 9) . Since the main contribution to the integral 
comes from from q ~ q c , one can replace W k — > v F q c s — 
2h c . Introducing u — uj — Wk, we have 



V2 



8tt N F v F \/^fK 
V2 q, 




-1/2 



(34) 



where u min = max{0,w - v F q c + 2h c } = 0, u max = 
min{ti;,ti; + v F q c + 2h c } = u> (we assume uj <C h c ). This 
expression takes a more transparent form if compared 
to the energy scale associated with the superconducting 
phase transition: 



Ac 



1.57 



Ao 
e F 



Au" 



1/2 



(35) 



r(fe,w) = -Im£f(fe,cj) 



We see that at the quantum phase transition into the 
LOFF state the Fermi-liquid behavior is destroyed by the 
superconducting fluctuations, and the magnitude of the 
fluctuation contribution to the decay rate is determined 
by the factor (A /e F ) 2 . 

In the generic case, when the inverse fluctuation prop- 
Im£^(<7, uj — Wk), (33) agator has minima at isolated points qc in the momen- 
tum space, see Eq. (|2U1) . the quasiparticle decay rate can 
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be written as the sum of the independent contributions 
from each minimum: 



r(fc,w) = 5^r (fc,«), 



(36) 



where 



r„ = 



N F 



d D i 



Here the integration is restricted by the condition < 
W k < u, and D — 3 or 2 is the dimensionality of the 
system. 

In three dimensions, after changing the variables q — 
qc — P, choosing the z axis in the momentum space 
along v k , and introducing u = uj — w a — \vk\p z , where 



"'a 



Vkq 



(«) 



2h c , one obtains: 

7 



4rt 2 N F \v h \ 



dp±PA 



du u 



K 2 [p\ + (u - uj + w a ) 2 /v 2 } 2 + 7 2 



In the limit uj — > 0, one can neglect u — uj in the first term 
in the denominator and calculate the integral over u: 



1 



87r 2 JVy|« fc | J 



x In 



dp±PA 



J 2 UJ 2 



K 2 (jP x -r >< „ • 



w 2 /v 2 ) 2 



(37) 



The result of the integration here essentially depends on 
whether w a is zero or not. 

If k is such that w a 0, then one can expand the loga- 
rithm in Eq. (|37() at uj — > and calculate the momentum 
integral. Substituting the result in Eq. (|36|l we obtain 



T(k,uj) 



1 



7\Vk\ 



tUJ 



(38) 



16tt 2 N f K 2 w 2 

w " — ^2 a w~ 2 . The fluctuation contribution to 
the quasiparticle decay rate has the energy dependence 
characteristic of the Fermi liquid, and its magnitude can 
be estimated as follows: 



r(fe,^) 

A 



Ao 



(39) 



On the other hand, if, for some k, one of u; a 's is zero, 
then w~ 2 diverges, making the expression (|38|l inappli- 
cable. Setting w a = in Eq. (|3Tjl . one obtains that the 
decay rate at uj — > is dominated by the contribution 
from the ath minimum: 



T(k,uj) 



1 



1 



16n N F \v k \K 



(40) 



so that 



Thus, the following picture emerges: In the generic case, 
the energy dependence of the quasiparticle decay rate 
due to the interaction with superconducting fluctuations 
turns out to be strongly anisotropic on the Fermi sur- 
face. While r oc uj 2 is almost everywhere, one has a 
non-Fermi-liquid behavior, r oc uj, on the lines defined 
by the intersection of the surfaces v^qi^ = 2h c with the 
Fermi surface. For the consistency of our calculation we 
have to assume that these intersection lines do exist, oth- 
erwise the coefficient 7 would be zero, see the discussion 
at the beginning of Sec. Ill CI 

Similar conclusions can be obtained in the generic 2D 
case, in which the decay rate has the same form 13611 . 
where, instead of Eq. 1)37(1 . we now have 



8ir 2 N F j\v k \ J_ 



dpi 



x In 



K 2 (pl+w 2 a /v 2 ) 2 



If w a 0, then 



T(k,uj) 



1 1< , 



(42) 



16tt N f K 2 w 



7,3 



(43) 



-7,-3 _ 



J2 a \ w a\ 3 - Using 7 - l/h c and K 



where 

v F /h 2 , we obtain 

T(k,uj) ^ Ao 
A ~ e F \A Q 

At the points on the Fermi surface where w a 
have the expression 



r(fc,a;) 



V2 



1 



47T N F \v k \y/jK 
whose magnitude can be estimated as 



1/2 



Ao 



Ao 



UJ 



1/2 



(44) 



0, we 



(45) 



(46) 



Similar to the generic 3D case, the fluctuation contri- 
bution to the decay rate strongly depends on the direc- 
tion of the Fermi momentum, showing non-Fermi-liquid 
behavior along some directions. Although the overall 
magnitude of the correction is larger than in the generic 
3D case, it is still proportional to the small parameter 
Ao/ep- It is interesting to note that the same frequency 
dependence of the decay rate can be obtained in the 
model of a nearly antiferromagnetic Fermi liquid in 2D, 
where the quasiparticle interaction with spin fluctuations 
becomes anomalously strong near some points, the "hot 
spots" , on the Fermi lineS 

Let us now compare our results with the decay rate 
at the second-order phase transition into the uniform su- 
perconducting state. Setting q c = and assuming an 
isotropic band dispersion, we have 
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instead of Eq. I|36|l . In the 3D case, repeating the calcu- 
lation steps leading to Eq. I|38[). one obtains: 



1 -yv F 2 

— Ul . 



64tt 2 N F K 2 hi 
Similarly, in the 2D case, 

1 jv F 



r( w ) 



128tt N F K 2 hl 



(48) 



(49) 



Thus, the Fermi-liquid character of quasiparticle excita- 
tions is not destroyed by the quantum fluctuations at the 
normal metal-to-uniform superconductor transition. 



IV. CONCLUSIONS AND OUTLOOK 

We have studied the order parameter fluctuations near 
the quantum phase transition at H = H c from the normal 
state to the LOFF superconducting state. We derived 
the general form of the fluctuation propagator C(q, v m ) 
at finite q and v m . In the systems suggested as good 
candidates for the experimental realization of the LOFF 
state, disorder is small or absent altogether. In the ab- 
sence of impurity effects, we analyzed the momentum 
and frequency dependence of the fluctuation propagator 
in both isotropic 3D and 2D cases, as well as in the case 
of generic spectrum. 

The fluctuation effects are more pronounced in the 
isotropic 3D case compared to the generic situation. This 
is because in the isotropic case the LOFF states are in- 
finitely degenerate leading to the large phase volume of 
fluctuations. The fluctuation contribution to the spin 
susceptibility diverges at H — > H c , 8\ <x (H — H c )~ x / 2 , 
and the quasiparticle decay rate shows a non- Fermi-liquid 
behavior, T(oj) oc uj 1 / 2 , at the quantum critical point. 
The magnitude of the fluctuation corrections is deter- 
mined by the parameter (Ao/ej?) 2 . Although this ratio 
is very small in conventional bulk superconductors, it can 
vary in a wide range in the atomic Fermi gases. It should 
be noted that our results rely on the assumption that the 
LOFF transition at T = is of second order. As men- 
tioned in the Introduction, this might not be the case in 



clean isotropic systems in the weak-coupling limit. How- 
ever, as the ratio Ao/eF grows so do the strong-coupling 
corrections to the mean-field free energy, which could sta- 
bilize the second-order LOFF transition. More theoreti- 
cal work is needed to check if this possibility can indeed 
be realized. 

In the generic case, which is expected to be applica- 
ble to crystalline paramagnetically-limited superconduc- 
tors, or to the atomic Fermi gases in optical lattices, the 
equilibrium wave vectors of the LOFF state form a set 
of isolated points in the momentum space. The phase 
volume of fluctuations is reduced, resulting in a non- 
singular spin susceptibility in the 3D case, and a weak, 
logarithmic, divergence of 5\ in the 2D case. Interest- 
ingly, the fluctuations in the generic case are still strong 
enough to cause the breakdown of the Fermi liquid at the 
quantum critical point, which manifests itself in a highly 
anisotropic energy dependence of the quasiparticle decay 
rate on the Fermi surface: T(oS) oc to on some lines in 3D, 
and F(o;) oc uj 1 ^ 2 at some points in 2D. We expect that 
the Fermi-liquid behavior, T(lu) oc lo 2 , will be restored 
throughout the Fermi surface away from the quantum 
critical point or in the presence of disorder, leaving the 
details for a future publication. 

There is a number of open questions concerning the 
assumptions we made and the effects we neglected in the 
present study, the order of the quantum phase transition 
into the LOFF state being particularly important. Even 
if the mean-field transition is of second order, this might 
no longer be the case if one takes into account the fluctu- 
ation renormalization of the free energy. Additional com- 
plications arise in realistic crystalline superconductors, in 
which the orbital effects and disorder should be included. 
Finally, it would be interesting to extend our calculation 
of the fluctuation corrections to nonzero temperatures in 
the critical region around the LOFF transition. 
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